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Abstract. A contribution of the electron-phonon interaction into the energy of a 
unidirectional charge ordered state (stripe phase) of two-dimensional electrons in GaAs 
heterostructures is analyzed. The dependence of the energy on the direction of the 
electron density modulation is calculated. It is shown that in electrons layers situated 
close to the (001) surface the interference between the piezoelectric and the deformation 
potential interaction causes a preferential orientation of the stripes along the [110] axis. 

PACS numbers: 73.20.Dx 
1. Introduction 

The observation of a resistance anisotropy of two-dimensional (2D) electrons at high 
Landau level occupancy ||, ||, |J is considered as an evidence for a formation of 
unidirectional charge ordered states (stripe phases) in such systems. These states 
were predicted in Refs. j5], |6], 0]. The Hartree-Fock calculations g |, 0] show that 
electrons at a topmost half-filled Landau level are separated in stripes of full and empty 
occupancy with a period of about several times the cyclotron radius. If the stripes are 
preferentially oriented along a certain direction, such states are expected to demonstrate 
the anisotropy in the longitudinal resistance: the low resistance along the stripes and 
the high resistance at the perpendicular direction ||, |J. 

The effect is observed in 2D electron layers in GaAs/AlGaAs heterostructures grown 
on (OOl)-oriented GaAs substrates. For such an experimental setup in the magnetic field 
perpendicular to the electron layer the low and high resistance directions are always 
pinned along certain crystallography axes of the host matrix, namely, [110] and [110], 
corresp ondingly. 

A mechanism that determines a preferential orientation of the stripes in GaAs 
heterostructures is not completely understood yet. It was suggested by Takhtamirov and 



Volkov [ID], that an anisotropy of the effective mass of 2D electrons may be responsible 
for the orientational pinning of the stripes. In the model developed in Ref. [|IIJ the 
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effective mass anisotropy is caused by an asymmetry of the quantum well potential 
confining the electrons to the (001) plane. This idea has been checked experimentally 
by Cooper et al [Tl\. In experiments |l f the influence of the symmetry of the quantum 



well potential on the resistance anisotropy has been studied. It has been found, that the 
resistance does not depend on the precise form of the confining potential, in particular, 
the same form of the anisotropy is observed for symmetrically confining 2D electrons. 
The authors of Ref. [|ll|] conclude that the asymmetry of the confining potential cannot 



be considered as an important factor. 

Another possible mechanism for the stripe orientation has been studies by the 
present author fl2]| . It was shown, that the piezoelectric interaction lowers the energy 
of the charge density wave aligned along either the [110] or [110] axes. Nevertheless, 
the piezoelectric mechanism alone does not explain why only one of two preferential 
orientations is realized. In this paper we address this problem. We extend the model 
proposed in |12| and take into account, together with the piezoelectric interaction, the 



deformation potential interaction. We find that, for the electron layers situated near the 
(001) surface of the sample, an interference between two channels of the electron-phonon 
coupling plays an important role in the stripe orientation. The mechanism considered 
explains the resistance anisotropy observed in the experiments. 

2. The model 

Let us consider a semi-infinite piezoelectric crystal that occupies a volume z < and 
contains an electron layer at the distance d from the surface. The static energy of the 
system can be presented in the form 

u= [ d 3 r( E£ + ^ )+ / sX + u^, (1) 



z<o 8tt 2 Jz>o 8tt 
where 

A = eEi - 4ir(3 iM Uki (2) 
is the electric displacement field, 

0~ik = KklmUlm + Pl,ik^l , (3) 

the stress tensor, E, the electric field, Uik, the strain tensor, \ikim, the elastic moduli 
tensor, fii^u the piezoelectric moduli tensor, e, the dielectric constant. To be more 
specific, we restrict our consideration to the case of a cubic lattice. The last term in ([]]) 
is the deformation potential interaction. It is chosen in the form 



Udef = J d rhp{u xx + u yy )5(z + d) , (4) 

where A is the deformation potential constant, p, the 2D electron density. Since we 

consider the model of an electron layer of zero-thickness, the interaction with u zz 
deformations is not included in (|J). 

The electric and elastic fields in Eq. ([[]) satisfy the following equations: 

VD = , (5) 
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(at z < 0), and 

VE = 

(at z > 0). 

At the free surface [z = 0) the boundary conditions are 
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At z = —d the normal component of the electric displacement field is discontinuous 
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The deformation potential interaction induces a tangential force applied to the medium 
in the z = —d plane. Rewriting the energy (HI) in the form 



(12) 



i=x,y 



where u is the displacement vector, we obtain the following expression for the force, 
applied to the unit area: 



(13) 



In equilibrium, this force is compensated by the stresses, and at z — —d the stress tensor 
satisfies the boundary condition 



(Ti 



z=-d+0 



= —Ft . (14) 

z= -d-0 

Using equations (0,1,0) and the boundary conditions one can reduce the energy (|l|) 
to the form 



U = - J d 2 r(epif - UiFi 



(15) 



where the electrostatic potential tp and the displacement field u are taken at z — —d. 
Their values are found from the solution of Eqs. (0-0) at the boundary conditions 
specified. 
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3. The stripe state energy anisotropy 



Let us calculate the energy fll^D for the stripe phase. We consider the system with 
the electron density modulated in a certain direction r s . The electron density can be 
presented clS cL Fourier series 

PM = E^ G ^PG„ , (16) 

Gn 

where G n = nq (n is integer), q || r s , \q\ = 2n/l, and / is the period of the stripe 
structure. For simplicity, we analyze the case of a unimodal charge density wave 

p(r P i) = po cos qr pl . (17) 



We calculate the energy (|l^) as a series in powers of the electron-phonon interaction 
constants: 

*L = U + U 2 + ... , (18) 
where S is the of the electron layer, 

the Coulomb energy in the absence the electron-phonon interaction, and the term U2 is 
quadratic in the interaction constants. For concrete physical systems considered below 
the contribution of the higher order terms in (|l^) is very small and it can be neglected. 
The term £7 2 determines the dependence of the energy on the direction of q. 

Let us first consider an isotropic crystal, for which the sound velocities do not 
depend on the direction of the sound propagation. For a cubic lattice this condition 
is realized if the elastic constants satisfy the relation en — C\2 — 2c^ = 0. For such a 
special case an analytical expression for U2 can be presented in a simple form. For the 
(001) electron layer the calculation of U2 gives the following result: 

U 2 = A + B cos 4ip + C sin 2ip , (20) 

where ip is the angle between q and the [100] axis. The angle dependence ( p0|) is 
determined by the parameters r\ = Cu/c^ and £ = qd. Using the strong inequality 
e ^> 1, that takes place in GaAs, we find for B and C the following approximate 
expressions 



B = E p 



, _ n_ _^ ( (, - »)W -y (*> - 3 >1 + t .[f + -{(,-!)])] , (21) 



c = e, e-^( '»" ° A y + emi - 1) - ^ - 3]) , (22) 

v 77 — 1 ' 



where E p = 9Tr 2 e 2 ple1 4 /64:e 2 Cnq, = 7r|e|/^ei4A/8ecn. We do not present here the 
expression for the parameter A which does not influence the orientation. 

In Eq. ( |2ll| ) the second term describes the anisotropy determined by the piezoelectric 
interaction, and the third term - the anisotropy caused by the interference of the 
piezoelectric and the deformation potential interaction. One can see that at £ — > 00 the 
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Figure 1. The parameter of the stripe energy anisotropy (Eq.(f22[)) in units of Ei for 
the isotropic crystal; r\ = 2.7, 2.5, 2.3 (from top to bottom). 

interference term tends to zero and the energy (p0|) obeys the C± v symmetry. If £ ~ 1 
the interference term is essential and the C± v symmetry is reduced to the C2 V one. We 
find that at rj < 2.7 the parameter B is positive for all £ and the global minimum is 
reached at ip m = 7r/4 or tp m = 37r/4 depending on the sign of the parameter C. At 
C < ip m = 71-/4 and the stripes are preferentially oriented along the [110] axis, while 
at C > ip m = 37r/4 and the [110] oriented stripe phase has the lowest energy. The 
dependence of C on d/l = £/27r is shown in Fig. [I]. One can see that at d/l > 0.3 the 
[110] orientation of the stripes is realized 

The results obtained are sensitive to the parameters of the system. Therefore, to 
determine the stripe orientation in GaAs heterostructures it is necessary to take into 
account the anisotropy of the elastic moduli. For this case we solve Eqs. with 
the boundary conditions specified, numerically. The following parameters are used for 
the calculations: c u = 12.3 • 10 10 N/m 2 , c 12 = 5.4 • 10 10 N/m 2 , c 44 = 6.0 ■ 10 10 N/m 2 , 
ei4 = 0.15 C/m 2 A = 7.4 eV, e = 12.5. The dependence of JJ% on the direction of q 
at I = 2 • 10 3 A and several d/l is shown in Fig. 0. In difference with the isotropic 
crystal, the minima of the energy are reached at q deviated from the [110] or [110] 
axes to the angles Atp ~ ±7r/12. If d/l is in the interval [0.23,1], the configurations 
with a small deviation of q from the [110] axis have the lowest energy. In this case 
one can expect that a domain structure is formed, and, in average, the system should 
demonstrate the minimum resistance in the [110] direction and the maximum resistance 
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Figure 2. The phonon contribution to the energy of the stripe structure (in 
units of E p ) for GaAs versus the direction of the electron density modulation at 
d/l = 1.0, 0.6, 0.5, 0.4, 0.25,0.23, 0.2 (from top to bottom). 

in the perpendicular direction. The calculations predict the largest resistance anisotropy 
at d/l = 0.4. At larger and at smaller d/l the anisotropy becomes weaker. At d/l > 1 
and d/l « 0.23 it disappears completely. At d/l < 0.23 the resistance anisotropy resets, 
but the high and the low resistance directions alternates. 

4. Discussion 

We have shown that the classical energy of a charge density wave in a 2D electron system 
embedded in a piezoelectric matrix depends on the direction of the wave vector. The 
effect is caused by the electron-phonon interaction. The minimum energy is reached at 
two or four different directions of the density modulation. Therefore, the stripes may 
form a polydomain structure. In the bulk samples polydomain structures may not show 
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any resistance anisotropy, while in 2D layers situated near the surface of the sample the 
resistance anisotropy should take place. 

We calculate the static energy of the stripe structure. Our results can be also 
understood as an effect of a virtual exchange of acoustic phonons between electrons. 
The electron-electron interaction in bulk isotropic piezoelectrics caused by the virtual 
exchange of acoustic phonons has been studied by Rashba and Sherman [I3~|. Our 



approach |L2| (see also Ref. ||14|1 , where the orientation of bi-layer Wigner crystals has 



been studied) reproduces the results of Ref. fll3| . Within such an interpretation the 
effect described in the present paper is a consequence of the virtual exchange by surface 
phonon modes. For the surface acoustic waves on the (001) surface the piezoelectric and 
deformation potential interactions give in-phase contributions to the matrix elements 
of the electron-phonon interaction jl5| , and the interference between two channel of 
electron-phonon interaction takes place. 

The unimodal approximation predicts that at d/l > 0.23 the monodomain stripe 
structure has the lowest energy for the stripes deviated from the [110] axis to the 
angle Aip « +7r/12 or Ai/j « —tt/12. In this case, the polydomain structure should 
demonstrate an anisotropic resistance with a minimum at the [110] direction. The 
largest resistance anisotropy is reached at d/l « 0.4. At d/l > 1 the anisotropy becomes 
exponentially weak. At d/l < 0.23 the anisotropy changes its sign. Note, that the last 
conclusion is specific for the unimodal approximation ( |T7j ) only. The higher harmonics 
in (|T6| ) will shift the transition point to a smaller value of d/l. 

We consider our model describes the orientation of stripe structures formed at high 
Landau levels. We analyze the anisotropy of the direct interaction between electrons. 
One can expect that the anisotropy of the exchange interaction is small, and it does not 
influence significantly the orientation. 

In the Hartree-Fock approximation the period of the electron density modulation I 
is approximately equal to 6£h, where £h is the magnetic length. In experiments |ll| the 
resistance anisotropy was observed at the magnetic field H « 2 T and for d m 2 • 10 3 
A. Therefore, for the Hartree-Fork I we find the ratio d/l ~ 2. At such d/l the surface 
effects are not important and only an exponentially small violation of the C± v symmetry 
may take place. Since in experiments this violation is quite large we suppose that the 
period I is large then it follows from the Hartree-Fock theory. 

It is interesting to evaluate the absolute value of the native anisotropy caused by the 
phonon mechanism. Using the parameters given before for d/l = 0.4, the filling factor 
v = 9/2, the electron density n = 2 ■ 10 11 cm -2 , and po ~ P (where p is the average 
density at the valence Landau level) we find the anisotropy energy E a « 0.7 mK per 
electron (we determine E a as the energy at q || [110] minus the energy at q || [110]). 

In our consideration we neglect the screening of the electron-phonon coupling caused 
by the polarization of the remote Landau levels. To evaluate the effect of screening in the 
unimodal approximation one can use the effective dielectric constant e e ff(q) = e(l+K q ) 
and the screened deformation potential A scr (q) = A/ (1+K q ), given by the random phase 
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approximation. Here 




,1 finer) ( $ 
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' n 2 P 2 



(23) 



In Eq. ( [23] ) f na = nF(e na ) is the Fermi factor, uj c is the cyclotron frequency, is the 
generalized Laguerre polynomial, and the prime on the sum excludes the valence Landau 
level (compare with Ref. |L6|]). Under such a substitution, the form of the dependence 
U2(ip) (Fig. |2|) remains almost unchanged, but the absolute value of the anisotropy is 
reduced as E^ rc « E a /(1 + K q ) 2 . Evaluation of the formula (p3|) for the parameters 
given above yields E S J C « 0.4 mK. 

The phonon contribution to the native anisotropy is comparable with the one given 
by the effective mass anisotropy mechanism [RJ (it gives the value of order of 1 mK). 
The experimental data for the anisotropy energy have been obtained in Ref. from 



the measurements of the resistivity in a tilted magnetic field. The data presented in [11 



are based on theoretical calculations by Jungwirth et al |16| . According to Ref. [jn| the 
anisotropy energy for the sample with a conventional heterostructure is higher (2.4 mK), 
then for the sample with a symmetric quantum well (0.5 mK). Thus, we conclude that 
two mechanisms of the anisotropy work in parallel, and the survival of the anisotropy 
in samples with symmetric quantum wells can be accounted for the phonon mechanism. 

The mechanism of the resistance anisotropy considered is this paper is essentially 
dependent on the distance between the surface of the sample and the electron layer. 
Therefore, it is desirable to investigate this dependence experimentally. Such a study 
may answer the question whether or not the surface effects play an important role in 
the stripe orientation. 

It is also of interest to investigate experimentally the influence of the electron layer 
orientation on the resistance anisotropy. To illustrate this point, we outline the results 
obtained for the (111) layer in the isotropic crystal. For such a system the energy (^C|) 
is modified to 

u 2 = A' + B'cos6tp (24) 

where ip is counted from the [Oil] axis. The coefficient B' is given by the following 
expression: 

10 / p" 2? 
B' = E p -(l- V - ^—-^[A( V 2 - v -l)(l + 20 + e(v + l)(fy - 7) 

+le(v+i)(v-s)-e(v-i) 2 }) (25) 

It is important to note, that the interference term in (|24]) does not depend on ip (it 
is included in A'). We find that the parameter B' is negative for all £, and rj > 2. 
Therefore, the minimum of the energy is reached at ipm = mr/3 (n is integer). It means, 
that monodomain stripe structures should demonstrate the low resistance along any of 
the [211], [121], or [112] direction, while the poly domain structures may not show any 
resistance anisotropy at all. Since this conclusion is not sensitive to the parameters 
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of the systems, the lattice anisotropy is not very essential in this case, and the same 

behavior is expected for the GaAs system. 
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